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Abstract 

m . 

Some multidimensional generalizations of the Fokker-Planck equation used 
by R. Friedrich and J. Peinke for the description of a turbulent cascade as a 
stochastic process of Markovian type, are considered. The exact solutions of the 
Cauchy problems for these equations are found with the operator methods. 

00 

^ 1 Introduction 

The understanding of the turbulence is one of the main unsolved problems of classi- 
cal physics, in spite of the more than 250 years of strong investigations initiated by 
D.Bernoulli and L.Euler. 

In the stochastic approach to turbulence the turbulent cascade is considered 
as a stochastic process, described by the probability distribution P(X, v), where A and v 
are the appropriate scaled length and the velocity increment respectively. Recently || 
R.Friedrich and J.Peinke presented experimental evidence that the probability density 
function P(X,v) obeys a Fokker-Planck equation (FPE) fU (see fig.l and fig.2 in |§): 



X 



dP{X,v) 
dX 
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.^ {XlV) + ^X,v) 



P(X,v), 



where the drift and duffusion coefficients D 1 and D 2 respectively are derived by analysis 
of experimental data of a fluid dynamical experiment (see fig.3 in [[|). 

In their paper Friedrich and Peinke use the following approximations for the drift 
and diffusion terms respectively 

D 1 = -av, a > 0; D 2 = cv 2 , c>0, 
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where a and c do not depend on A. 

In our previous paper ||, using the method of M.Suzuki || and the Feynman's 
disentangling techniques we gave the exact solution of the Cauchy problem for the 
Eq. (|ip with more realistic approximations (see fig. 3 in ||) for D 1 and D 2 : 

D 1 = -a{X)v, a(A)>0; D 2 = c{X)v 2 , c(A) > 0. (2) 

In Section 2 of this paper we consider the Cauchy problem for the following n- 
dimensional generalization of the Eqs. ([I]), (0): 

BP 

— = b (X)P(X, v) + &!(A)v • V V P(A, v) + c(A)v 2 A v P(A, v), (3) 

P(0,v) = p(v). 
This equation is derived from the n-dimensional FPE 

§ = -V v • D '(A, v)P(A, v) + ± ^-D U K v)P(A, v) (4) 

with the drift vector D 1 (A,v) = — a(A)v and the diffusion tensor P 2 (A, v) = c(A)v 2 l. 
( Consequently 6 (A) = n[a(X) + 2c(A)j and 61(A) = a(A) + 4c(A) ). 
In Section 3 we find the exact solution of the Cauchy problem 

BP 

— = 6 (A)P(A, v) + 6i(A)v ■ V V P(A, v) + c(A)(v • V V ) 2 P(A, v), (5) 

P(0,v) = p(v). 

The Eq. @ is derived from the FPE (|) with the drift vector D X (A, v) = — a(A)v and 
the diffusion tensor D 2 (X, v) = c(A)vv. ( Consequently 60(A) = na(X) + {n 2 + n)c(A) 
and 61(A) = a(A) + (2n + l)c(A) ). 

Because of the analytic expressions of the diffusion tensor, one may regard the cases 
(||) an d (H) respectively as the "isotropic" and "degenerate anisotropic" (det D = 0) 
ones. 

2 Exact Solution of the Cauchy Problem (g) 



In the spirit of the operational methods ( see ||— |iT| ) we have for the solution of the 
problem (|3|) 

P(A, v) = (^exp + £ [60(a) + 6 a (s)v ■ V v + c(s)v 2 A v ] ds^j tp(y), (6) 
where the symbol exp + J x C(s)ds designates the V.Volterra ordered exponential 
exp + [ X C(s)ds = l+\imJ2 [ X dAi dX 2 . . . dX k C(X 1 )C(X 2 ) . . . C(X k ). (7) 

Jo n ~* oo ^"'0 Jo Jo 
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The linearity of the integral and the explicit form of the operators in (|6|) permit to 
write the solution P(X, v) in terms of the usual, not ordered, operator valued exponent 



P(A,v) =e WA; e 



where for convenience we have denoted 



A) (A) „/3i(A)v.V v +7(A)v 2 A v 



tp(v) 



ft(A)=/ 6,(a)da, (j=0,l); 7 (A) = / c(s)ds. 



(9) 



Consequently (from now on ""' means 



&(0) = 0, (3'JX) — bj(X), (J = 0,1); 7 (0) =0, 7 '(A) = c(A). (10) 



Since 

A(A)v-V v , 7 (A)v 2 A v 
the solution (M) can be written in the form: 



P(A,v) =e tt,w e 



^o(A) o/ 3 1 (A)v-V Vo7 (A)v 2 A, 



,/3o(A) o7 (A)v 2 A Vo /3 1 (A)vV v 



V(v) (12) 



To write the expression ([12] ) in a final form we will use the following formulae 
for acting of the pseudodifferential operators [|T^| - e /3v Vv and e TV Av on arbitrary 
functions of v: 

e /?(A)v.v v/ ( v ) = f ( w e m) (13) 



and (for r(A) > 0) 



e r(A)v 2 A v ^ (v) ; (> 



r(A)A e 



/OO s 2 
Q-^mg f ve (n-2)r(X)±s\ ds 
-oo ^ ' 



/47rr(A) 

Here Ag is the operator of Laplace-Beltrami on the sphere Si in lZ n ||15||— (IT 

A, 



(14) 



d 2 n — 1 d 1 { 



<9i> 



A^ = A^, \ = -l(l + n-2), 1 = 0,1,2,..., 



(15) 



(16) 



A; — 1,2, 3, — (ij^, <i; jn 



(2/ + n-2)(n + /-3) 



(n-2)!Z! 

The eigenfunctions Y^(Q) are the spherical (harmonic) functions on Si which consti- 
tute the orthonormal basis in L 2 (Si), i.e for any function g(v) we have: 



s(v)^M)=EEa^>) 

1=0 k=l 
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9iM = {^\g)= [Y^(9)g(v,e)dS u 



(17) 



Si 



where 



v = {\v\,6 1 ,6 2 ---6 n „ 1 ) = {v,6), 6 ir --e n - 2 E [0,7r], n _ x 6 [0, 2tt] , (18) 

v± = v cos 9i , 
t>2 = f sin 9\ cos 

^n-l = v S111 #1 " - " srn $n-2 COS Q n -\ , 

v n = v sin Q\ • • • sin 6 n _ 2 sin n _i , 

c!5i = sin"" 2 9 1 sin"" 1 2 • ■ • sin #„_ 2 d#id0 2 ■ ■ • d0 n _ x . 
Thus for the exact solution of the problem (§) we obtain from the Eg . (|T2|) 

e /3o(A) 



P(A,v) 



>/4tt7(A) 



^(A) 



OO a 2 

e JtPI 

-oo 



,7(A)A e 



OO d'l.n 



•^T^(ve^W + ( n - 2 ^ ±s )d S 



(19) 



(A)+(n-2) 7 (A)± 



Z=0 fc=l 



') 



ds, 



where v e 1Z n , /3 (A), /?i(A) and 7(A) are from (§), and A/ and c?; n are from 
In partucular, for v G TZ 3 , we have (6*i = 9 and #2 = <p) 

{ 9<p ~ sin 9 89 Sm 89 + sin 2 9 dip 2 ) 



P(A,v) 



p /3o(A) . ^00 Q 2 

( ,7(A)A e / e -3^)^ (ve ft(A)+7{A)±s ) ds 

-00 ^ ' 



,/9o(A) 



OO s 2 
— OO 



V / 4^(A) 

00 1 



(20) 



i=0 m=-Z 



ds, 



where Y{ n {9, 93) are the Laplace functions on the 1-sphere in P 3 . 

Substituting the expression (|i~9"|) in the Eq. (|3|) one can see immediately that P(A, v) 
is a solution of the problem (||) and, according to the Cauchy theorem, it is the only 
classical solution of this problem. 

3 Exact solution of the problem (|5|) 

Proceeding as in Section 2, we have for the solution of the problem (||): 



P(X, V) = e /3o(A) e /3i(A)vVv+7(A)(v-V v )' 



(21) 
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where /3 (X), /3 2 (A) and 7(A) are from (|]). 
Since 

"ft(A)v-V v , 7 (A) (v- V v ) 2 ] =0, (22) 
we can factorize in the Eq. (f2l|): 

p(A,v) = e^e^^e^^Vv) = e^e^^V^^v). (23) 
Now using (|13"D and the formula 

e 7W(vV v ) 2 ^ v ) = 1 Z 00 e -47W^ (ve ±s ) ds, (24) 

we find for the exact solution of the problem (^): 

/3 (A) .qo s 2 

P(A, v) = -== / e"^)^ ( ve ftW±») dSj (25) 

where /3o(A), /3i (A) and 7(A) are defined in (H). 

Substituting the expression ([25|) in the Eq.(|5]) one can check that P(X, v) is solution 
of the problem (||), and, according to the Cauchy theorem, it is the only classical 
solution of the problem. 

4 Conclusions 

The exact solutions of the Cauchy problems (^) and are obtained, using the 
disentangling techniques of R. Feynman and M. Suzuki's methods for solving the 
1-dimensional Fokker-Planck equation. The problems (|3|) and ([]) are two of many 
n-dimensional generalizations of the 1-dimensional FPE, used by Friedrich and Peinke 
in their description of a turbulent cascade as a stochastic process of marcovian type. 
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